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AF0$R-77-33O3  FINAL  TECHNICAL  REPORT 
1  APRIL,  1977  TO  31  JANUARY.  1982 

1 )  Abstract 

We  have  provided  the  only  significant  experimental  and  analytical 
data  on  pseudo-stationary  oblique-shock-wave  reflections  in  monatomic 
(Ar),  diatomic  (O2,  N2,  air),  tri atomic  (CO2)  and  polytomic  (SFg) 
gases  in  the  form  of  isopycnics  and  shock  Mach  number  -  wedge  angle 
(M  ,0  )  -  plots  of  value  to  the  military  (AFWL,  DNA,  ARO,  BRL.NSWL, 

3  W 

ORES,  etc.),  industry  (R  &  D  Associates,  Physics  International,  NASA, 
etc.)  and  universities.  Our  laboratory  results  have  been  invaluable 
to  field  experiments  and  to  computational  fluid  dynamicists.  The 
latter  have  had  to  develop  new,  improved  codes  to  match  our  accurate 
experimental  data. 

Our  results  on  the  production  of  neutrons  and  y-rays  by  explosive- 
driven  implosions  in  D2  have  now  been  published.  They  have  attracted 
an  international  interest.  The  method  provides  a  simple  means  of 
studying  fusion  plasmas  at  extreme  temperatures  and  pressures  as  well 
as  solid-phase  transitions  from  carbon  to  diamond  and  other  new  materials. 

The  reports  on  the  dynamics  of  dusty-gas  flows  produced  by  shock 
waves  are  "best-sellers"  and  are  out  of  print.  We  hope  to  test  our 
analyses  in  our  new  10cm  x  20cm  dusty  gas  shock  tube  built  for  this 
purpose  during  the  year  with  funds  from  the  Canadian  Defence  Research 
Establishment  Suffield  (ORES).  The  analytical  and  experimental  data 
has  important  applications  to  combustion,  dust  explosions,  damage  from 
blast  waves  and  cosmic  gasdynamics. 

Our  extensions  of  the  random  choice  method  (RCM)  have  aided  us 
(and  others)  greatly  in  analyzing  explosion-implosion  dynamics, 
detonations,  dusty-gas  flows,  and  viscous  heat-conducting  vibrational ly- 
excited  shock-wave  transitions  in  air  for  weak  spherical  N-waves. 


The  virtue  of  this  method  is  that  it  does  not  suffer  from  artificial 
viscosity  (implicit  or  explicit)  smearing  of  shock  waves  and  contact 
surfaces.  If  some  method  could  be  found  to  apply  RCfl  to  pseudo-  and 
non-stationary  flows,  it  would  revolutionize  computational  shock  fluid 
dynamics. 

Our  analytical  v/ork  on  swirling  turbulent  combustion  flows  has 
progressed  very  well.  Experimental  data  are  now  being  accumulated 
using  laser  Doppler  anemometry  (LDA).  This  work  is  not  only  of 
interest  to  researchers  in  various  establishments  (including  AFOSR) 
but  to  industry  as  well.  In  this  regard,  we  have  received  some 
financial  support  from  Canadian  Pratt  and  Whitney  to  set  up  the  LDA 
equipment. 

2 )  Research  Objectives 

The  research  objectives  are  to  make  significant  advances  in  the 
field  of  "Aerophysical  Aspects  of  Gas  and  Plasma  Flows’’,  useful  to 
the  varying  supporting  agencies  and  to  research  and  development  in 
industry.  Our  objectives  have  been  met  to  a  large  extent.  This  state¬ 
ment  can  be  verified  from  our  numerous  publications  in  distinguished 
journals  and  by  our  peers  in  the  field  internationally. 

3 )  Detailed  Summary  and  Perspectives 

Although  the  phenomenon  of  oblique-shock-wave  reflections  has  been 
known  for  a  long  time,  it  is  only  within  the  period  of  this  Grant  that 
order  and  understanding  has  been  brought  to  this  subject.  Once  and  for 
all,  we  have  shown  that  the  domains  and  boundaries  of  regular  (RR),  single 
Mach  (SMR),  complex-Mach  (CUR)  and  double-Mach  reflections  (DMR)  can  be 
plotted  as  a  function  of  the  initial  conditions,  namely,  shock  wave  Mach 
number  M5  and  wedge  angle  Oy,  -plane,  for  a  perfect  gas.  If  real-gas 
effects  take  place  then  the  initial  temperature  Tq  and  pressure  Po  must 
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be  added.  The  results  for  argon  (monatomic  gas)  are  shown  in  Fig.  1. 

It  is  seen  that  the  analysis  is  good  and  agrees  with  experiment  for  a 
perfect  gas  in  the  range  tested  for  Ms<10.  The  results  for  air,  N2  and 
O2  (diatomic  gases)  appear  in  Fig.  2.  It  can  be  seen  that  the  agreement 
with  a  perfect  gas  in  this  case  is  quite  good  except  for  the  SHRJCMR 
transition  line.  The  experimental  data  consists  of  our  own  work  and  from 
others.  Figure  3  shows  the  results  for  CO2  (a  triatomic  linear  molecule). 
Again  the  data  agree  well  with  analysis  in  the  range  l<Ms<10. 

These  two-dimensional  data  along  with  interferometric  plots  of 
lines  of  constant  density  (isopycnics  have  formed  the  basis  for  testing 
the  validity  of  several  computational  methods  of  predicting  the  flow 
quantities.  So  far,  the  computational  schemes  are  quite  good  for 
predicting  shock  shapes  and  wall -density  distributions  but  poor  for 
giving  results  on  isopycnic  shapes  and  distributions.  These  disagreements 
are  being  investigated  by  many  government,  university  and  industrial 
establishments  in  a  number  of  countries.  Undoubtedly  better  computational 
procedures  will  evolve  in  the  near  future. 

The  foregoing  laboratory  results  have  been  successfully  applied 
by  AFWL,  DNA  and  BRL  among  others  to  predict  shock-wave  configurations 
and  pressure  loadings  in  field  tests  on  vehicles  and  missile  sites. 

The  transition  lines  themselves  must  be  fine-tuned  in  order  to 
improve  the  agreement  between  analysis  and  laboratory  experiments.  A 
great  deal  has  already  been  done  in  this  respect  and  this  will  be  the 
subject  of  future  reports  under  the  new  Grant  AF-AFOSR-82-0096. 

The  production  of  fusion  plasmas  in  deuterium  cannot  be  described 
more  simply  than  in  the  brief  paper  in  Physics  of  Fluids  by  Glass  and 
Sagie,  "Application  of  Explosive-Driven  Implosions  to  Fusion",  (Ref.  1). 

A  copy  is  enclosed.  This  unique  method  of  producing  neutrons,  yrays 
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and  solid-phase  transitions  (graphite  to  diamond  of  10-20um  dia) 
offers  simple,  new  possibilities  for  many  avenues  of  physical  research 
of  interest  to  industry,  universities  and  the  military.  As  a  matter 
of  fact,  3M  of  Canada,  who  are  interested  in  industrial  diamonds,  have 
taken  up  our  research  and  development  program  on  the  use  of  explosive- 
driven  implosions  to  generate  diamonds  from  graphite.  Several  R  &  D 
personnel  from  3M  are  stationed  at  our  Institute  pursuing  this  v/ork. 

It  is  an  excellent  example  of  technology  transfer  from  a  university 
to  industry. 

Although  the  dynamics  of  dusty-gas  flows  is  of  much  importance 
to  science,  industry  and  the  military,  very  limited  research  and  develop¬ 
ment  was  conducted  in  the  past  in  this  field.  In  a  more  practical 
vein,  very  little  is  known  about  nonstationary  drag  and  heat  transfer 
experienced  by  individual  particles  accelerated  by  shock  or  blast  waves. 
Neither  is  it  known  precisely  i/hat  impact  pressure  such  flows  impose 
on  structures,  nor  what  damage  may  result  to  vehicles,  hardened  missile 
sites  or  buildings  affected  by  high-speed  high-pressure  dusty-gas  flows. 

In  order  to  understand  such  problems,  the  two  avenues  of  analysis  and 
experiment  have  been  undertaken.  Our  analytical  work  on  the  flow  in 
a  dusty-gas  shock  tube  (Figs.  4  to  6)  and  the  passage  of  a  shock  wave 
through  a  dusty-gas  layer  (Figs.  7  and  8)  (in  great  demand  internationally), 
will  form  the  basis  for  our  experimental  work.  For  the  latter,  we  have 
built  and  will  instrument  a  3-3/4  x  8  in.  dusty-gas  shock  tube  60  ft. 
long,  in  order  to  check  our  analyses  by  using  glass  spheres  as  homogeneous 
dust  in  the  lO-BOum  range.  Optical  methods,  laser  Doppler  anemometry, 
pressure  gauges,  impact  gauges,  heat-transfer  gauges  and  other  devices 
will  be  used  for  this  purpose.  We  estimate  the  dusty-gas  shock  tube 
will  cost  $500,000.00  and  it  is  being  supported  by  Canadian  funds.  Here 
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is  an  excellent  example  of  AFOSR  benefits  coming  from  an  almost 
entirely  Canadian  funded  project. 

The  random-choice  method  (RCM)  of  analyzing  nonstationary  planar, 
cylindrical  and  spherical  flows  with  shock  waves  has  proven  to  be 
one  of  the  most  pov/erful  and  accurate  computational  methods  available 
today.  Its  power  and  accuracy  lies  in  the  fact  that  it  makes  no  use 
of  finite  difference  methods  which  suffer  from  an  inherent  explicit  or 
implicit  artificial  viscosity  that  smears  shock  fronts  and  contact 
surfaces  in  a  flow.  Instead  it  solves  the  Riernann  problem  at  each  time 
step  in  the  numerical  analysis  in  a  random  fashion.  The  results  provide 
sharp  fronted  shocks  and  contact  surfaces.  This  method  has  been  applied 
successfully  at  our  Institute  to  numerous  problems  such  as  explosion, 
implosion  and  detonation-wave  dynamics,  nonlinear  planar-wave  interactions 
and  dusty-gas  dynamics.  However,  we  have  also  succeeded  in  using  the 
method  for  spherical  N  .vave  shock-front  transitions  with  viscosity, 
heat  Conductivity  and  vibrational  excitation.  The  N-waves  were 
generated  by  exploding  wires.  The  explosion  process  was  modelled  by 
considering  the  H-waves  produced  by  a  small  pressurized  sphere 
which  is  suddenly  ruptured.  The  sphere  d'smeters  and  the  initial 
conditions  determine  the  properties  of  the  N-waves  that  simulate  those 
actually  obtained  from  exploding  wires.  Some  of  the  results  appear  in 
Figs.  9  and  10.  The  comparison  between  analysis  and  experiment  is 
good.  This  research  will  help  us  understand  the  properties  of  N-waves 
generated  in  the  atmosphere  by  SST's  and  will  explain  the  order-of- 
magnitude  discrepancy  between  measured  and  predicted  risetimes  of  sonic 
boom  N-waves.  The  risetime  controls  the  human-startle  effect.  The 
shorter  the  risetime  (microseconds),  the  more  annoying  is  the  boom. 

So  far  the  RCM  has  been  successful  in  treating  nonstationary  one- 
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dimensional  flows.  What  is  required  for  the  treatment  of  pseudo-stationary 
two-dimensional  oblique-shock-wave  reflections  or  nonstationary  spherical 
shock-wave  reflections  is  a  RCM  algoritliiJi  to  handle  such  complex  flows. 

So  far,  despite  the  endeavours  of  many  able  computational  fluid 
dynami cists,  the  search  has  proved  to  be  very  elusive. 

The  analytical-numerical  program  with  coupled  gasdynamic  and 
chemical -dynamic  equations  of  motion  to  deal  with  turbulent,  svnrling, 
combusting  flows  has  been  completed.  Now  emphasis  has  been  placed  on 
laser  Doppler  anemometer  (LDA)  experimentation  shakedowns  in  order 
to  measure  the  turbulence  components  and  Reynolds  stresses  at  a  given 
station  of  a  combustor.  The  results  can  then  be  used  as  data  input  in 
the  analysis  for  predicting  these  quantities  at  any  other  station 
downstream.  Their  new  measurements  can  be  made  at  the  predicted  stations 
in  order  to  check  the  various  models  of  turbulence  and  the  coupled 
equations  of  motion. 

As  a  fall  out  of  the  LDA  work,  it  has  been  confirmed  that  particle 
size  and  their  velocity  can  readily  be  measured  to  lib.  This  will  make  it 
possible  to  measure  fuel-spray  size  and  velocity  distributions  in 
combustors  which  is  of  much  importance  to  the  jet-engine  industry  in 
the  USA  and  Canada. 

Finally,  our  research  and  development  work  in  the  foregoing  areas 
is  continuing  in  a  vigorous  manner.  Many  reports  are  in  preparation 
in  the  five  major  areas  considered.  Several  papers  have  been  sent  to 
distinguished  journals.  One  has  already  been  accepted  for  publication 
in  the  Proceedings  of  the  Royal  Society  of  London  and  others  will  follow. 
Details  will  be  given  under  Grant  AF-AFOSR-82-0096  in  the  near  future. 

Dr.  I .  I .  Glass 
Principal  Investigator 
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FIG.  1  EXPERIMENTAL  VERIFICATION  OF  OBLIQUE  SHOCK-WAVE-REFLECTION 

REGIONS  AND  BOUNDARIES  IN  THE  (Ms,Ow)-PLANE.  PERFECT  MONATOMIC 
GAS  Y  =  1.667,  LAW  (1970),  CMR  IN  HELIUM  (?);  LAW  AND  GLASS 
(1971)  CMR  IN  ARGON  (V),  BAZHENOVA  ET  AL.  (1976)  DMR  IN  ARGON  (O). 
PRESENT  RESULTS  IN  ARGON:  D,  DMR;  T,  CMR;  A,  SMR;  #,  RR. 
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FIG.  2  EXPERIMENTAL  VERIFICATION  OF  OBLIQUE  SHOCK-WAVE-REFLECTION 

REGIONS  AND  BOUNDARIES  IN  THE  (Ms,Ov;)-PLANE.  PERFECT  DIATOMIC 

GAS  Y  =  1.400.  - REAL  GAS.  t,  C,  AIR  (DATA  FROM  SMITH 

1945);  V,  9,  AIR  (WHITE  1951);  O  ,  V,  A,  0,  OXYGEN  (LAW 
AND  GLASS  1971);  B  ,  V,  AIR  (BAZtIENOVA  ET  AL.  1976);  A, 
NITROGEN  (BAZHENOVA  ET  AL.  1976);a  ,  T,  A.  0,  NITROGEN  (BEN-DOR 
AND  GLASS  1978). 
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Fig. 6  Variation  of  shock  v/ave  Mach  number  Mg  with  diaphragm 

pressure  ratio  Pl,  .  - - -  Ms  ^^sed  on  frozen  speed  of 

sound  ap^, - Mg  based  on  equilibrium  speed  of  sound 
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FIGURE  8.  Flow  quantities  at  t  =  3.  12xl0~^  s  (  $).  (a)  Pressure, 

(b)  Mass  concentration,  (c)  Temperature,  (cl)  Velocity.  -  gas, 

particles,  (i)  Beginning  of  formation  of  a  reflected  shock  wave, 
(ii)  Beginning  of  formation  of  a  contact  region. 
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Application  of  explosive-driven  implosions  to  fusion 

1. 1.  Glass  and  D.  Sagie“’ 

lor  A^mpace  Siudiei.  L  nuvnny  of  Toronto,  Toronto,  Canada 
iRcccived  H  May  I'JSl,  accepted  20  November  I'JSl) 

Our  ecpiosive-driven  implocion  facility  was  used  to  prtxluce  hemispherical  implosions  in  a  stoichiometric 
mixture  of  deutenum-oxygen  A  high-resolution  scintillator  detection  system  measured  neutrons  and  y  rays 
resiiilitig  mr>sl  likely  from  the  fusion  of  deuterium. 


I.  INTRODUCTION 

The  Institute  for  Aerospace  Studies  hemispherical  im¬ 
plosion  chamber  is  a  unique  device  for  producing  ex¬ 
plosive-driven  implosions  which  are  stable  and  well- 
focused  at  the  geometric  center  in  a  safe  and  reusable 
facility.'  ft  has  been  used  as  a  driver  for  launching  hy- 
pervelocity  projectiles,^  to  generate  intense  planar 
shock  wave.s,  *  and  to  produce  diamonds  from  graphite.' 
In  the  present  .study  is  was  utilized  to  bring  about  fusion 
in  a  deuterium- oxygen  plasma. 


II.  EXPERIMENTAL  EQUIPMENT 

The  implosion  chamber  (Fig.  1)  consists  essentially  of 
two  massive  mating  steel  jilates.  The  rear  iiluie  con¬ 
tains  a  20-cm  diar.i  hemispherical  cavity.  Into  it  i.s  lu¬ 
ted  an  explosive  package  consisting  of  a  copper  liner  to 
which  IS  bonded  a  shell  of  suiierline  PETN  secondary 
explosive  (of  about  3  mm  in  thickness  weighing  97  g  and 
releasing  about  0.6  MJ  of  energy  for  the  case  reported). 
The  front  plate  contains  an  exploding  nickel  wire  (0.13- 
mm  diamxi-mm  long)  and  the  gas  inlet  and  outlet. 

Both  plates  are  fastened  together  by  32  bolts.  The  hem¬ 
ispherical  cavity  is  filled  with  a  stoichiometric  mixture 
of  deuterium-oxygen  (in  this  case  about  55  atm  releas¬ 
ing  about  0.4  MJ  of  energy).  The  gas  is  detonated  at  the 
geometric  center  by  the  exploding  wire.  The  gaseous 
detonation  wave  instantly  and  simultaneously  explodes 
the  PETN  on  impact,  thereby  generating  a  well-focused 
implosion  wave.  This  wave  reflects  at  the  geometric 
center,  leaving  a  small  pocket  of  plasma  at  extreme 
pressure  (niegabars)  and  temperature  (millions  of  de¬ 
grees)  in  which  a  deuterium  fusion  reaction  occurs.  To 
our  knowledge  this  is  the  most  direct  method  of  initiat¬ 
ing  fusion  using  only  chemical  energy.  We  have  tried 
other  indirect  methods  to  obtain  fusion,  where  a  small 
capsule  containing  deuterium  was  placed  at  the  focal 
point  of  the  implosion.  One  of  these  was  successful  and 
provided  almost  identical  results,  thereby  lending  sup¬ 
port  to  the  simple  direct  technique.  Details  can  be 
found  in  Ref.  5. 

The  detection  system  for  sensing  neutrons  and  y  rays 
(produced  by  the  neutrons  and  their  interactions  with  the 
steel  implosion  chamber)  consists  of  two  scintillator- 
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photoiiiultijilier  assemblies  (Fig.  1).  The  lirsl  detector 
IS  located  at  the  outer  surface  of  the  front  plate,  30  cm 
Iroiii  the  implosion  focus  while  the  second  is  BO  cm 
from  the  focus.  The  first  oscilloscope  displays  the  en¬ 
tire  ignition-detonalion-implosion  process  lasting  about 
50  trsec  and  the  subsequent  events.  The  second  oscillo- 
scojie  IS  designed  to  display  the  full  undisturbed  shape 
of  the  first  signal  coming  from  'he  second  detector  and 
therefore  sweeps  at  a  rate  of  50  '■<'  :/division.  From 
Its  shape  it  would  have  been  possible  to  obtain  the  neu¬ 
tron-velocity  distribution  and  its  flux.  However,  the 
second  oscilloscope  can  only  trigger  if  a  large  enough 
signal  (-0.5  V)  is  produced  by  the  first  detector.  This 
requires  about  two  neutrons  to  cross  the  scintillator 
within  10  nsec.  The  threshold  level  is  essential  to  jire- 
vent  false  triggering  arising  Irum  the  photomulli|)lier 
dark  curreni,  cosmic  ravs  and  ignition  noise. 

III.  EXPERIMENTAL  RESULTS 

Records  of  voltage  versus  tune  from  the  first  detector 
lor  two  runs  without  ia)  and  with  fusion  (b)  are  shown  in 
Fig.  2.  Initially,  there  are  large  oscillations  arising 
from  the  capacitor  discharge  to  the  exploding  wire, 
which  are  damped  out  in  about  35  nsec.  In  the  case  of 
no  fusion  when  a  stoichiometric  mixture  of  Hj-Oj  is  de¬ 
tonated,  no  other  signals  appear.  However,  with  a 
stoichiometric  mixture  of  Dj-O,  when  fusion  occurs  at 
about  50  psec,  about  20  negative  signals  appear  in  a 
random  lime  and  amplitude  distribution  over  a  period  of 
about  50  psec.  The  maximum  amplitude  of  about  0.2  V 
corresponds  to  a  single  impact  by  a  neutron  or  a  > 
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FIG.  1.  Schemafic  of  implosion-rhaTnbcr  facility  and  scintil¬ 
lator  detectors. 
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I-'ICi,  2,  ()Sfill()Si'ope  rpc'orr)  from  doloctor  No.  1,  (a)  With¬ 
out  fusion,  211, ■■  0_,  mixture  tit  27.2  atm  and  127  g  PETN  ex¬ 
plosive  (hi  W'ith  fusion,  2n;i-  O,  mixture  at  u-t.-t  atm  and  !>7 
i;  I’ETN  explosive;  li)  arrival  of  implosion  and  beginning  of 
events,  lii)  ignition  noise. 


photon  in  the  MeV  range,  obtained  from  a  calibration 
of  the  scintillators  at  the  University  of  Rochester.  In 
this  case  the  second  oscilloscope  did  not  trigger  as 
none  of  the  events  generated  a  large  enough  signal.  We 
estimate  that  the  total  yield  was  a  few  thousand  neu¬ 
trons. 

The  distribution  of  events  can  be  explained  by  the 
scattering  of  neutrons  as  they  encountered  the  steel 
chamber  walls.  Iron  has  a  total  cross  section  of  about 
3.5  b  at  2.45  MeV.'’  Other  components  of  the  steel  alloy 
like  Mn  and  Co  also  contribute  to  the  neutron  scattering. 
Inelastic  scattering  by  iron  atoms  produces  mainly  0.85 
MeV  and  1.24  MeV  y  rays.  However,  neutron  capture 
by  Fe'"  produces  metastable  Fe’’ with} -ray  energies  up 
to  10  MeV.  The  total  number  of  neutrons  and  photons 
generated  per  reacting  neutron  is  about  three.  The  cou¬ 
pling  angle  between  the  chamber  and  the  detector  is  a- 
bout  tenfold  greater  for  the  scattered  neutrons  and  > 
rays  than  for  the  neutrons  generated  at  the  implosion 
focus.  Consequently,  the  scintillator  is  far  more  effi¬ 
cient  in  detecting  the  indirect  scattered  neutrons  and  > 


rays  than  the  direct  neutrons  from  the  implosion  focus. 
Owing  to  the  large  attenuation  of  the  inelastically  scat¬ 
tered  neutrons,  the  phenomenon  is  spread  out  and  de¬ 
layed  as  recorded  in  Fig.  2(b). 

IV.  DISCUSSION  AND  CONCLUSIONS 

A  consideration  of  the  plasma  parameters  shows  that 
the  peak  temperatures  may  be  obtained  at  a  radius  of  a- 
bout  10  pm  with  an  ion  density  of  about  5x10^^  ions/cnA 
The  implosion  time  is  then  about  10''“  sec.  Although 
radiation  and  conduction  heat  losses  are  significant, 
they  would  allow  temperatures  to  be  reached  up  to  a  few 
keV. 

Discussions  of  different  work,  with  greater  complex¬ 
ity,  to  produce  neutrons  in  deuterium  by  explosive 
means  can  be  found  in  Ref.  5,  as  well  as  some  consid¬ 
erations  of  the  possibilities  of  scaling  the  present  ap¬ 
paratus  to  obtain  thermonuclear  fusion.  Here,  we  have 
shown  that  neutrons  and  y  rays  can  be  obtained  from 
nuclear  reactions  by  very  direct  means  from  an  explo¬ 
sive-driven  hemispherical-implosion  focus  in  D2-O2 
mixtures.  There  is  little  doubt  that  temperatures  in 
the  keV  range  were  reached  thereby  approaching 
thermonuclear  fusion  conditions.  Much  work  remains 
to  determine  the  details  of  the  physical  properties  of 
such  plasmas  and  the  resulting  nuclear  collision  pro¬ 
cesses. 
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